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Abstract 

An algebraic model based on Lie-algebraic techniques is applied to the analysis of 
thermodynamic vibrational properties of diatomic molecules. The local anharmonic 
effects are described by a Morse-like potential and corresponding anharmonic bosons 
associated with the U (2) algebra. A vibrational high temperature partition function 
and the related thermodynamic functions are derived in terms of the parameters 
of the model. The thermal expansion operator is defined and obtained using the 
renormalized frequency of the model. 

1. Anharmonic Algebraic Model 

Algebraic models have been used very successfully in nuclear physics and have led to new 
insights into the nature of complex many body systems [1]. Similar methods have been 
proposed to describe molecular excitations [2,3]. The methods combine Lie algebraic 
techniques, describing the interatomic interactions, with discrete symmetry techniques 
associated with the molecules. In the framework of the algebraic model [2], the anhar- 
monic effects of the local interactions are described by substituting the local harmonic 
potentials by Morse-like potential, associated with the U (2) algebra. The one-dimensional 
Morse Hamiltonian is written in terms of the generators of Z7(2), 
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where A is a constant. The eigenstates, | [N],v], correspond to the U(2) D SU(2) 
symmetry-adapted basis, where N is the total number of bosons fixed by the potential 
shape, and v is the number of quanta in the oscillator, v = 1, 2, . . . , [^] . 
The anharmonic effects are described by anharmonic boson operators [2], 



b = 



- Jz 
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where v is the Morse phonon operator with an eigenvalue v. The operators satisfy the 
commutation relations, 
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The harmonic limit is obtained when N — > oo, in which case 
the usual boson commutation relations. 
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The Morse Hamiltonian can be written in terms of the operators b and b 
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which corresponds to vibrational energies 
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where uj is the harmonic oscillator frequency. The Morse potential spectrum thus leads 
to a deformation of the harmonic oscillator algebra. A more detailed relationship between 
the Morse coordinates and momenta and the SU (2) generators can be derived through 
a comparison of their matrix elements [5]. Note that for an infinite potential depth, 
N —> oo, the Morse potential cannot be distinguished from the harmonic potential. 

The anharmonic model has been developed to analyze molecular vibrational spectra 
[2-8]. It provides a systematic procedure for studying vibrational excitations in a simple 
form by describing the stretching and bending modes in a unified scheme based on U(2) 
algebras which incorporate the anharmonicity at the local level. 

The aim of this paper is to apply the algebraic approach to the vibrational high- 
temperature thermodynamics of a diatomic molecule and obtain the basic thermody- 
namic functions in terms of the parameters of the anharmonic model. These results must 
be combined with the translational and rotational thermodynamic functions in order to 
compare with experiment, in a fashion similar to that of reference [9]. Here, we con- 
centrate our attention on the derivation of analytic results which will be applied in a 
subsequent paper [10]. 



2. Vibrational Partition Function 

The vibrational partition function of a diatomic anharmonic molecule is 
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[N/2] 

Z N =J2 ^ (6) 

v=0 

where j3 — l/fe^T, the vibrational energies e v are given in equation (5) and N is the fixed 
total number of anharmonic bosons. 

Introducing new parameters, a = l3h ^ J " , N = [^] and I = [^] — v, the vibrational 
partition function can be written as, 

N 

Z N = e- a J2e-^ N2 °- l2 \ (7) 

1=0 

At high temperatures T, for N large and a small, the sum can be replaced by the 
integral, 




Z N = J ^e-^o+i) I e s* ds 



/aNo 

J e s2 ds (8) 



where s = yj^l- This integral can be evaluated exactly in terms of the error function, 
crfi(VoTJVo) (sec [11]), 



7 1 N TT 



^e-^+^erf^v^iVo). (9) 



Equation (9) represents the high-temperature value of the vibrational partition func- 
tion in the Morse- like spectrum. When No — > oo, aNo ^ 1> the harmonic limit of the 
model is obtained, 



Zoo ~ ~ = e 2k s T (10) 

which coincides with the harmonic vibrational partition function of a diatomic molecule 
at high temperatures. 

Having the partition function Zn in terms of the parameters of the algebraic model, 
we are now in position to derive the basic thermodynamic functions. An algebraic ap- 
proach has been used in [12] to study the thermodynamic properties of molecules. How- 
ever, the partition function in [12] uses an approximation of the classical density of states, 
while we have derived an explicit function in terms of the parameters of the algebraic 
model which can now be used to compute the high-temperature thermodynamic functions 
in analytic form. 
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3. Thermodynamic Vibrational Functions 

3.1. Mean Vibrational Energy 

The mean vibrational energy is given by 



Um = — - — In Zm = — . (11) 

df3 2Z N da 1 ; 



Taking into account that 

^~£-<* + .>*» + ^ 

we obtain the following expression for the mean vibrational energy in terms of the parti- 
tion function Zn, 



Substituting Zm by expression (9) gives, 



2 y u/ 2 



(14) 



The harmonic limit is obtained from equation (13), when iVo — > oo and Z^v is given 
by (10), 

tfoo ~ ^ + M (15) 
This is the classical mean energy of a diatomic molecule at very high temperatures. 

3.2. Specific Heat 

The vibrational part of the specific heat is, 

8U N Tiujq dU N 

Substituting Un with equation (13) and using (12), we obtain 



k B k B N e- a ( 1 N e- , 

CN = - + -2Z N -{ N ° a -2-^Z N -) (17) 

This represents the vibrational specific heat in the anharmonic model in terms of the 
partition function Zn- All anharmonic contributions depend on the temperature. The 
explicit value is obtained from equation (17) by substituting Zn by expression (9). When 
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TVo — ► oo, the harmonic limit of the model gives the vibrational specific heat of a diatomic 
molecule at very high temperatures, 

Oc ~ k B . (18) 

3.3. Mean Number of Anharmonic Bosons 

The mean vibrational energy in the anharmonic model can be written in terms of mean 
number (vn) of anharmonic quanta, each with energy hu>o, 

U N = Hujo (Vat) + (19) 
Substituting Un by equation (13), we obtain (i/jv) in terms of the partition function Zjy, 

N I N e- a 

{vn) = - + ^-j^- (20) 

Using expression (9) in equation (20), we obtain the high-temperature value, 
N Q k B T [Nq~ fk B T\? e N °^ 

<-> = - + 2^-VvUrJ e^/gy- (21) 

The harmonic limit is obtained from equation (20) when Nq — ► oo and Zn is given 
by expression (10), 

<"»> ~ < 22 » 

3.4. Free Energy 

The free vibrational energy in terms of the partition function Zn is given by 

F N = ~\nZ N (23) 

Substituting Zn with equation (9) gives the free vibrational energy in the algebraic 
model at high temperatures, 



Using expression (10) in equation (23), we obtain the classical harmonic result for the 
free vibrational energy at very high temperatures, 

F^^ksT^. (25) 
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4. Thermal Expansion 

It is well-known that the thermal expansion in a molecule is a result of the local an- 
harmonic effects and cannot be explained by the harmonic models. It depends on the 
positions of the atoms and their displacement from the equilibrium. As an approxima- 
tion, we can write the coordinate operator x in terms of the anharmonic operators b and 
6 f (see also [5]) as 



*-te( t+,t ) (26) 



Here, u> depends on the number of quanta in the oscillator [4-7] , 



N N 

Now, we define the thermal expansion operator 



. N + l 2v\ 
^ = ^Kr-77 (27) 



^=V^( S+St+ ^) (28) 

where represents the displacement from the equilibrium position. The harmonic limit 
is obtained when N —> oo, in which case u> — > ujq an d —> 0. 

Evaluating the mean value of the operator, Xn, between the states | [N]v) and taking 
into account (27), we obtain the mean thermal expansion Xm in terms of the parameters 
of the anharmonic model, 



Xn = \I^- ^n(n + i) i'-^h) (29) 

The anharmonic effects are essential at high temperatures T. When N^>v, equation 
(29) becomes, 



2/iu>o N 

Taking into account the high-temperature value of the mean number of quanta (22), 



X»~J^l(l + P^). (31) 



2^uj N \ huj N 

The thermal expansion of anharmonic diatomic molecule has been obtained in [13], 
using a Morse-type potential and classical methods. The comparison of equations (29) 
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and (30)-(31) with the results in [13] at high temperatures shows very good agreement. It 
also shows that is proportional to the mean separation between the atoms at T = and 
that iV 2 is proportional to the depth of the potential well, the latter being in agreement 
with the algebraic model [4-7]. 

5. Conclusion 

We have applied the algebraic model to those thermodynamic properties of diatomic 
molecule which at high temperature strongly depend on its anharmonicity. We have de- 
rived the vibrational partition function and the related thermodynamic functions, such as 
mean vibrational energy and specific heat, in terms of the parameters of the model. The 
next step is to consider the rotational excitations of the molecule. The thermal expan- 
sion, which is a typical anharmonic effect, has been derived in terms of the parameters 
of the model and its high-temperature value is in a good agreement with classical anhar- 
monic results. The model can be further applied to the thermodynamics of polyatomic 
molecules. This and other applications will be discussed elsewhere. 
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